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ABSTRACT 



o 

^ ■ We present follow-up observations of the far-infrared (FIR) sources at 90, 150 and 

\ 180 iLtm detected as part of the ISOPHOT EBL project, which has recently measured 

the absolute surface brightness of the cosmic infrared background radiation (CIRB) 
for the first time independently from COBE data. We have observed the fields at the 
North Galactic Pole region in the optical and near-IR, and complement these data 
with SDSS photometry, and spectroscopy where available, and present identifications 
^Nf^ I of the 25 FIR sources which reach down to ~ 150 mJy in all three ISOPHOT bands. 

. Identifications are done by means of full spectral energy density fitting to all sources in 

■ the FIR error circle areas. Approximately 80 per cent are identified as star- forming or 

star-bursting galaxies at z < 0.3. We also find that more than half of the counterparts 
have disturbed morphologies, with signs of past or present interactions. However, only 

■ 20 per cent of all the sources are uniquely matched with a single galaxy - 40 per cent 
, are blends of two or more of these nearby star- forming galaxies, while another 20 per 

cent are likely blends of nearby and fainter galaxies. The final 20 per cent are likely 
to be more luminous IR galaxies at higher redshifts. The blended sources have an 
. effect on the FIR source counts. In particular, taking into account realistic confusion 

^ • or blending of sources, the differential FIR counts move down by a factor of ~ 1.5 and 

I steepen in the 100 to 400 mJy range. 

Key words: Galaxies: evolution - Infrared: galaxies - galaxies: starburst 
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1 INTRODUCTION 

The initial detection of the cosmic infrared background 
(CIRB) in the fa r-infrared (FIR) usi n g COBE bet^y e en 10 
and 240 |j.m by ' Hauser et all (llQQSlVlFixsen et alj (|l998h . 
and prehminarily by Puee t et al J(|l996), showed that a large 
part of all the radiation produced in the universe since the 
recombination epoch is contained in the CIRB. In the con- 
text of understanding galaxy formation and evolution, and 
indeed the history of all luminous matter in the universe, it 
is thus important to understand what and where the CIRB 
contributors are. This has been a subject of intense work 
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ever s ince the COBE results (for rev iews see lHauser &: Dwe^ 
I2OOII : iLagache. Puget k Dolell2005l ). 

The ISOPHOT EBL project aimed to determine the 
level of CIRB independently from the COBE data by making 
use of the better resolution of ISO which allows looking into 
the very deepest regions on the FIR sky in-between Galactic 
cirrus clou ds. The FIR source counts from this project were 
reported in ljuvela. Mattila &: Lemkd (|2000l ). The CIRB sur- 
face brightness in three wavelengths 90, 150 and 180 ^im was 
recently presented in Juvela et al. (2009) - this constituted 
the first independent verification of the absolute CIRB level 
since the various COBE results. In between 150 and 180 ^im 
the CIRB was measured to be 1.1 =b 0.3 MJy sr~^ and at 90 
^im a 2(7 upper limit of 2.3 MJy sr~^ was measured, all val- 
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ues consistent with the COBE measurements obtained with 
the DIRBE and FIRAS instruments. 

However, only the bright tail of FIR sources contribut- 
ing to the CIRB can be resolved with data from past mis- 
sions, reaching to typical levels of somewhat fainter than 
^ 100 mJy at 180 /xm. The FIR source counts at 90 to 
180 /im in our project con tribute ^ 10 — 20 per cent of t he 
CIRB surface brightness (|juvela, Mattila Le mke 2000). 
This result is similar to other FIR studies (e.g . ^Dole et al.., 
l200ll : lFraver et al.l[2006bl : iMatsuura et aLllioOTH when look- 
ing at the FIR sources directly using ISO^ Spitzer, and 
AKARL Statistically one can probe fainter FIR popula- 
tions by using e. g. stacking metho ds, as shown bv lDole et al.l 
lioOS) and Pauovich et al.l (|2007l l . Due to significantly less 
severe confusion, the resolved fraction of individual sources 
in the mid-IR 15 and 24 iJ,m mid-IR wavelengths is much 
highe r, around 60^to 75 per cent (e.g. Elbaz et al. 2002; 
IPapovich et al.1 l2004l : iRodighiero et al.1 l2006l ). iFraver et all 
(|2006al ) report resolving 60 per cent of the CIRB also at 
70 fxm using fluctuation studies and extrapolated source 
counts. At wavelengths longer than this, up until ^ 200 ^tm, 
the fraction of directly resolved CIRB sources has essentially 
not changed since the various ISOPHOT studies. The recent 
BLAST experiment reports resolving the FIR background at 
500 ^tm b y stacki ng a Spitzer 24 \im catalog against their 
FIR map (D evlin et al.l l2009). 

During the decade of FIR follow-up studies it has also 
become evident how difficult it is to identify th e optical and 
near-IR counterparts of these FIR sources (e.g. Pat r i s et al. 
2003 : Tavlor et al. 2005; Ovabu et al. 2005; Dennefel d et al. 
20051 : ISaiina et al.ll2006l : lFraver et al.ll2006bl) . Typically only 
a third, or a half, of the FIR sources have unambiguous coun- 
terparts at shorter wavelengths. The problem is this: because 
of the large beam sizes there are many possible counter- 
parts in optical images within the positional error circles 
of FIR sources. Hence, to aid identification other waveband 
data, especially MIR or radio data, and often quite elabo- 
rate decision trees are typi cally used to merge catalogs with 
uniqu e counterparts (e.g. Tavlor et al.l 1200 5': Ovab u et al.l 
Another difficulty arising from the large beam sizes 
is the high level of confusion among the FIR detections re- 
sulting in blended objects, which in turn make the optical 
identifications even more complicated. These difficulties in- 
directly affect mid-IR studies as well, since many interesting 
cosmological or galaxy evolution results from the MIR sur- 
veys require information of the spectral energy distributions 
(SEDs) of the FIR part of the spectra of the targets. These 
are just not available to similar depths, unless the MIR data 
itself is used to identify sources and construct FIR SEDs, 
leading into possibly unhealthy dependencies. 

Nevertheless, a basic expectation and consensus has 
been built where the FIR sources are either local IR-bright 
sources dominated by cold cirrus emission or higher red- 
shift sources whose strong IR-luminosity raises them into 
the class of luminous or ultra-luminous IR-galaxies (LIRGs 
or ULIRGs,> lO^^L© or > lO^^L©, respectively). The rela- 
tive amounts of sources in these classes depend on the depth 
of the survey (see e.g. Rowan- Robinson 2009). To interpret 
the cosmological evolution of IR-bright galaxies, it remains 
crucial to understand the exact nature of galaxies emitting 
at 100 to 200 ^im where the CIRB peaks. The imminent 
Herschel data will probe these galaxies to further distances 



(|Pilbrattll20Q8l ) and it is important to determine the char- 
acteristics of the more nearby bright FIR galaxy population 
for reference. 

In this paper we revisit the source detections of the 
ISOPHOT EBL project by presenting optical and NIR 
follow-up observations of the ISOPHOT NGP fields. Rather 
than trying to force individual optical and NIR counterparts 
for each FIR source, we assess carefully the possibility of 
confused, or blended, sources. In addition, we do not use 
mid-IR or radio data to identify sources. We present the 
properties of all the unambiguous matches and discuss the 
level of likely confusion in the FIR source list. If it were to 
be shown that a significant fraction of ISOPHOT sources 
are actually confused sums of several objects, or even whole 
galaxy groups or clusters, FIR source counts and their evo- 
lutionary interpretation might have to be revised. 

In a second paper of the series (Vaisanen et al., in prepa- 
ration) we will concentrate on the correlations of the FIR 
source positions with significant peaks in the galaxy surface 
density, any signs of clustering and any bona fide galaxy 
clusters. 

We use a spatially flat cosmology with Hq — 71 km s~^ 
Mpc~\ Q = 1, and Qm = 0.27 throughout. 



2 DATA 

2.1 Observations and data reductions 

The ISOPHOT observ a tions are fully described in 
ljuvela. Mattila Lemk^ (|2000l ). The North Galactic Pole 
(NGP) field with an area of 1.07 sq.deg covers 2/3 of the 
whole ISOPHOT EBL project area and contains 22 individ- 
ual sources in the 170 and 150 \im. bands, and 25 sources at 
90 /im. Sources were detected down to ^ 150 mJy in all three 
bands, while simulations indicated a typical 70 per cent or 
more completeness at ~ 200 mJy. The three 'extra' detec- 
tions in the 90 |am are in fact sources very close to another 90 
|j,m source, and the two of them share the same single longer 
wavelength ISOPHOT detection; we therefore have followup 
of 22 different fields-of-view. All these fields were targeted 
with ground based observations presented in this paper. Ob- 
servations were performed using the 2.56-m Nordic Optical 
Telescope (NOT) on La Palma. Optical data were taken 
on 8-10 April 2003 using the ALFOSC instrument, and the 
near-IR observations were carried out on 13-15 April 2003 
using the NOTCam instrument. 

ALFOSC is a 2x2k CCD camera, giving a field of view 
of 6.4'x6.4'. We obtained 900 sec R and /-band (filters 76 
and 12) images of the target fields, split into three separate 
offset integrations. The data were reduced and combined in 
the standard manner using IRAF. 

NOTCam is a Ixlk "Hawaii" near infrared camera. The 
low resolution mode of the instrument was used, giving a 
4'x4' field of view. We obtained iC-band images of all the 
fields. Target fields were observed in 1-minute dithers for a 
total of 25-minute integration time per field. The flatfielding 
(twilight flats were used), sky frame construction, sky sub- 
traction, and registering and co-adding of the frames into 
final NIR images was performed using the XDIMSUM pack- 
age along with dedicated IDL routines. 

Seeing remained in the range 0.8 — 1.0 arcsec in iT-band 
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and 0.9 — 1.3 in the optical bands throughout our observing 
runs. The astrometry for both optical and infrared images 
was obtained through matching GSC2 stars - we estimate 
the accuracy to be better than OA'' throughout the frames. 

We note that while the NGP area has not been ob- 
served by other targeted FIR surveys in the past, it will 
shortly be covered by the Herschel Thousand Degree Sur- 
vey (PLEalesQ) 

2.2 Photometry 

All phot ometry was derived usin g;; the SExtractor (v2.3.2) 
software (|Bertin &: Arnout'3l 19961 ). Both aperture and 'auto' 
magnitudes were used as detailed below. We maximised the 
detection efficiency by performing systematic and compre- 
hensive testing of the SExtractor settings with our data. 
Numbers of real and potentially spurious objects were de- 
termined by using both deeper data, i.e. checking the deeper 
optical images when performing NIR photometry, and neg- 
atives of the final images in case of optical frames. 

We calculated numerous aperture magnitudes, as well 
as the AUTO and BEST magnitudes, which dynamically 
adjust the size of the aperture depending on the shape and 
size of the object. By using curves- of- growth with bright 
non-saturated stars, we determined aperture corrections for 
slightly different seeing conditions and fixed size apertures. 
When giving colours we typically refer to magnitudes cal- 
culated inside 2.6" diameter apertures, approximately 2.5 
times typical seeing, corrected to the K-hand seeing value. 
However, total magnitudes, adopting the BEST magnitude 
which we find to be most robust in our uncrowded fields 
with both extended galaxies brighter than R—23 mag and 
compact objects throughout our magnitude range, are used 
when a single-band brightness is given for an object. We 
checked these to be totally consistent with the total magni- 
tudes derived from curves of growth used for standard stars. 
An exception is made with bright galaxies clearly larger than 
^ 5" for which these total magnitudes are used throughout, 
also in the SEDs. In summary, we estimate our photometry 
to be accurate to approximately 0.03 mags in the optical and 
0.05 in the NIR data. Finally, we checked that our photome- 
try was consistent within few percent with 2MASS galaxies 
found in our areas. 

Zero-point calibrations were obtained through observa- 
tions of standard stars throughout the nights, which were 
photometric. The absolute calibration to Vega-based system 
is estimated to be accurate to 0.02 and 0.05 mag in the op- 
tical and NIR, respectively. The extinction coefficents were 
fixed to average values of the nights, to 0.08, 0.05, and 0.10 
in /, and iT-bands, respectively. In addition, a colour- 
term of 0.1 X {R — I) was used to derive the final optical 
magnitudes. 

Since detection of optical objects was done with co- 
added R-\- 1 frames, the R and / band source lists have the 
same number of objects. These were cross-correlated with 
the K-hand source Ust to produce source lists of 3-band de- 
tections, optical detections only (since RI data is deeper 
than NIR), and NIR-only (including very red objects). 



2.3 SDSS data 

We also extracted uqr iz photometry from the SDSS DR7 
(|Abazaiian et al.ll2009l ) archives for all our fields to comple- 
ment our own dataset. SDSS data are shallower than our 
own optical RI data, but they extend the wavelength range 
of the brighter galaxies and thus make SED fitting signifi- 
cantly more reliable. After testing with different magnitudes 
we found that we acquire the best results, the smoothest 
combined SEDs, when adopting total magnitudes, as de- 
scribed above, for our own photometry, and Petrosian mag- 
nitudes for the SDSS galaxies brighter than r < 20 mag, and 
SDSS model magnitudes for those fainter than this limit0. 



3 ISOPHOT SOURCES IN NGP FIELDS 

3.1 Matching of FIR and optical/NIR sources 

As discussed in Section [1] identification of FIR sources is dif- 
ficult because of confusion and the large beam sizes. Typical 
approaches include using supplementary MIR or radio sur- 
veys in helping to decide which of the many galaxies within 
the beam is the correct counterpart, and then running like- 
lihood ratio methods taking into account the distance of 
counterpart candidates from the FIR source as well as the 
surface densities of all relevant populations. Often the one 
MIR and/or radio galaxy giving the largest likelihood is 
chosen as the counterpart, unless no candidat es are found 
above a given probabil i ty threshold (see e.g. iMann et al.l 
I2OO2I : lOvabu et aLlliooBI : iTavlor et al.ll2005l l. Further analy- 
sis, SED fitting etc., is based on this identification. While the 
advantage obviously is to have a well-defined and quantita- 
tive method, the downsides include possible biases depend- 
ing on the "indicator source" and especially the difficulty 
of assessing the possibility of multiple galaxies contributing 
to the FIR source. In this paper we attempt a much more 
straight-forward philosophy, though admittedly one less ro- 
bustly defined. 

We first searched for all sources with R, /, and i^-band 
detections within 60'' of the centre of the ISOPHOT de- 
tection. Making use of the SExtractor star vs. galaxy neural 
network classifier, we excluded clear stellar objects by select- 
ing only those with a mean CLASS parameter less than 0.75 
in all the three bands. Typically there remain 20 to 50 ex- 
tended objects per ISOPHOT error circle. The optical/NIR 
sources were then matched also with the SDSS catalogs by 
searching for common objects within l'' radius. 



3.1.1 SED fitting and other analysis of candidates 

Next we fitted template SE Ds using the HYPERZ soft- 
ware teolzonella et al ] I2OOOI ) to all sources found in the 
ISOPHOT error circles. The whole range of wavelengths 
from UV to FIR was fitted, our optical and NIR data- 
points along with the five SDSS bands and three ISOPHOT 
points. The fainter targets not having SDSS photometry 
were fitted with RIK-\-lSO only. The model templates were 
taken from, or calculated using, the GRASIL library/code 



1 http://herschel.esac.esa.int/Docs/KPOT/KPOT_accepted.html 

#KPOT_seales01_2 ^ see " htt p : / / www . sdss . org / D R7 / algorithms / photometry. html | 
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(jSilva et alJ Il998l ) - in particular we used a range of Sc- 
type SEDs of different age, different age ellipticals, fully 
evolved Sa and Sb galaxies, as well as model fits to the 
SEDs of NGC 6946, NGC 6090, M82, Arp 220, and HRIO, 
which respectively represent star- forming dwarfs, interact- 
ing/starbursting galaxies, warmer starbursts, (ultra) lumi- 
nous IR galaxies, and extremely obscured objects. 

When fitting the templates within HYPERZ, we turned 
off reddening since the GRASIL SED models already have 
dust and gas properties with consistent extinction and re- 
emission mechanisms and effects built in throughout the 
wide spectral range. We assumed a ^ 0.15 mag uncertainty 
in the relative matching of FIR and ground-based parts of 
the SED due to uncertain aperture effects. 

The source catalogs in the target fields, now associated 
with photometric redshifts, best-fit SED templates, and es- 
pecially the related chi-square values and probabilities at- 
tached to the HYPERZ fits, form the basis of the selection 
of counterparts. Note especially that we do not use any in- 
formation of the distance to the FIR detection centre in 
the analysis, but treat all sources in the error circle equally. 
[Rodighiero et al. (2005) cite a typical blended FIR source 
distance of 55^^ in their 90 ^tm data, and such cases would be 
severely misinterpreted by e.g. just picki ng the closest cand i- 
date to the FIR detection centre (see also lSaiina et al.ll2006l ). 
On the other hand, we note that by fitting known templates 
to objects we cannot get away from the bias against totally 
unknown classes of sources, which plagues all FIR identifi- 
cation studies. 

We also searched the SDSS spectroscopic archive for 
any available spectra for the bright galaxies in the target 
fields. Within the ISOPHOT error circles 12 redshifts were 
found, and for these objects the redshift was fixed in the 
HYPERZ fit. 

To have a quantitative measure of the morphology of the 
brighter galaxies in the field, we used the GALFIT program 
(|Peng et al.ll2002l l to model the surface brightness distribu- 
tions of the RI combined images. 

3.1.2 Breakdown of the candidates 

Based on the types of optical-NIR objects within the 1 ar- 
cmin search radius characterized by their magnitudes, SEDs, 
model fits and redshifts (photometric or spectroscopic), as 
well as morphologies, we separated the 22 FIR objects in the 
present sample into two groups: In the first instance, we find 
5 fields which have a single bright K ^ 14 mag galaxy in the 
ISOPHOT detection circle. All these galaxies are more than 
two magnitudes brighter than any other object in the area 
and they are all well fit with a template SED all the way 
from UV and optical to FIR. They also have SDSS spec- 
troscopy available which is consistent with their nature as 
FIR sources. They are all brighter than i^=14.8 mag and we 
note that using aver age bright K-hand galaxy counts (e.g. 
IVaisanen et ahlboOOl l. there is less than 10 per cent probabil- 
ity to find a K < 14.8 mag galaxy inside a 60^^ radius circle. 
We refer to these five unambiguous ISOPHOT counterparts 
as the "bright sample" . 

The second group of ISOPHOT targets, 17 of them, is 
called the "faint sample", referring to the brightness of the 
counterpart galaxies, not that of the FIR flux. There are 
essential differences amongst these sources, however. More 



than half of the fields do actually have fairly bright K ^ 
13 — 16 mag galaxies in them, but there are usually two to 
four of them, suggesting possible blending. The rest of the 
fields, on the other hand, look initially empty of any obvious 
counterpart candidates. When investigating the faint sample 
individually below, we categorize them according to how and 
why the cases are ambiguous. 



3.1.3 Further modelling of SEDs 

To have an independent characterization of the brighter 
counterparts with redshifts, or secure estimates thereof, 
we modelled them wit h th e radiative transfer models of 
Efstathiou et al. (2000") and Efstathi ou Rowan- RobinsonI 
(2003) (hereafter ERR03). The mo dels incorporate the stel - 
lar population synthe sis model of iBruzual Chariot! (|l993h 
and the dust model of lSiebenmorgen Kriigel (|l992h which 
treats small transiently heated grains and PAH molecules 
as well as large classical grains. The models have been 
shown to be in good agreement with the spectral en- 
ergy dlstribiry£ns_^ 

les jEfstathiou et al.1 l2000l; [ Efstathiou Rowan- RobinsonI 



120031 : iFarrah et al.ll2003l : iTavlor et alll2005l ) 



3.2 Unambiguous bright galaxies 

All the FIR sources of the bright sample defined above 
have FIR source quality fia^ values of ^ ^ 3 (see 
ljuvela. Mattila Lemk^ l2000h in at least one detected 
band. The sources are listed in Table [1] along with some 
characteristics based on GRASIL model fits, whereas Table[2] 
lists the ERR03 model fits. Based on our optical imaging, all 
these sources are disk galaxies. Images from the combined 
R + I ALFOSC imaging are shown in Fig. □ SEDs of all 
sources within 60^^ radius from the ISOPHOT source are 
shown in Fig. [2] and the best-fitting SED to the bright can- 
didate is also overplotted. The detailed description of each 
of the five fields, with discussion leading to the identification 
of the source, its SED and morphological fit, is presented in 
the Appendix. 

As can be seen from the Tables [1] and [2] the likely coun- 
terparts are generally nearby z < 0.2 (strongly) star- forming 
late type spiral galaxies. The general characteristics of the 
population is discussed below in Section [4. II 



3.3 Confused fields 

There are 17 cases in the present FIR sample where we can- 
not find a single unambigious counterpart, though in many 
cases there are fairly bright, K = 13 — 16 mag galaxies in the 
vicinity. In fact, in more than half of the fields there are two 
to five K < 16 mag galaxies within the ISOPHOT detection 
error circle, while statistically only 0.4 K < 16 galaxies are 
expected. It is likely that the majority of these cases are 
indeed confused FIR sources where the IR flux originates 
from several individual sources. Images of fields based on 
ALFOSC data are shown in Fig. [3] We discuss the cases in- 
dividually in the Appendix, while a summary is presented 
in Table [3l and the general physical characteristics are dis- 
cussed in Section 112] Figure g] shows the RIK-\-FlR SEDs 
of the target fields; note that the overplotted SED templates 
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NGP01 



NGP03 



NGP20 




Figure 1. Optical R -\- I combined images of the bright sample of the ISOPHOT NGP source fields. The large circle of one arcmin 
radius shows the location of the ISO detected source. The small red circle identifies the likely optical counterpart discussed in the text 
and zoomed-up in the inset - its best-fit SED template is shown in Fig. [2] North is up and East left. 

Table 1. Bright unambiguous ISOPHOT counterparts are listed with their SDSS redshifts, the best fitting UV-to-FIR GRASIL SED 
template from HYPERZ, the absolute K magnitude, and effective radius in kpc from GALFIT. L/^ is the log of IR luminosity, in 
solar units , integ rated over 8 — 1000 ia.m of the best-fit GRASIL model and SFR is calculated from this using the FIR/SFR relation of 
iKennicuttI (|l998l V 



FIR source 


counterpart 


z 


SED 


Mk 


LiR 


Ts 


SFR 


comments 


NGPOl 


2MASSXJ13412738+4042166 


0.088 


Sc" 


-24.22 


10.94 


2.9 


15 


disturbed 


NGP03 


2MASSXJ13422216+4022017 


0.131 


SB^ 


-24.72 


11.17 


6.3 


26 


interacting 


NGP20 


J134934.9+390730.1^ 


0.140 


Sc^ 


-24.28 


10.92 


2.9 


14 


undisturbed 


NGP22 


SDSSJ135054.71+385847.2 


0.086 


Sc" 


-23.42 


10.62 


3.6 


7 


disturbed 


NGP24 


UGC 08793 


0.0081 


Scd^ 


-19.63 


8.49 


2.2 


0.05 


warped, many bright HII regions 



^ nominally 15 Gyr Sc-model 
b NGC 6090 model 
^ uncatalogued - see text 
d NGC 6946 model 



are indicative of typical templates and typical redshifts, and 
not the best fits for the individual cases. 



4 DISCUSSION 



4.1 The nature of unambiguous counterparts 



We found in Section [3. 2 1 that five of our 22 FIR sources, i.e. 
23 per cent, are unambiguously identified with low-redshift 
z < 0.2 star-forming spiral galaxies. Apart from NGP24 all 
are just about at or above the LIRG luminosity criterion 
of > 10^^ Lq. Apart from NGP20 all have evidence of dis- 
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Table 2. The bright ISOPHOT counterparts fitted with models of lEfstathiou et al.l (|2000') : 'Efst athiou Rowan- RobinsonI (|2003l V The 
logs of the cirrus and starburst IR luminosity components, Lc and Lsb-, respectively, in solar units, do not necessarily equal the total 
listed in Table [T] based on GRASIL models, though the results are clearly consistent in each case. All starburst models assume an age of 
26 Myr and an initial Ay of the molecular clouds of 50. We also list the resulting peak star-formation rate, and the log of dust mass in 
solar units and the temperature in K. 



FIR source 


z 


Lsb 


Lc 




SFR 


dust 


log Kust 


NGPOl 


0.088 


10.66 


11.16 


1.3 


20 


17 


8.3 


NGP03 


0.131 


10.68 


11.15 


1.1 


22 


25 


7.5 


NGP20 


0.140 


10.62 


11.07 


1.3 


20 


17 


8.3 


NGP22 


0.086 


10.55 


10.95 


1.3 


17 


17 


8.2 


NGP24 


0.0081 




9.16 


0.3 




17 


5.9 



^ The temperature of the dominant 0.24yLim graphite grain species in the best-fitting model in the FIR is listed. 

^ Dust mass is estimated as in iTavlor et al.l (|2005l ) but the method has been generalised to take into account the distribution of dust 
temperature predicted by the model. 




1.0 10 100 1,0 10 100 1.0 10 100 

X (yL6m) 
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Figure 2. SEDs of bright individual galaxies in CIRB ISOPHOT fields. Circles represent our RIK photometry within 60^^of the FIR 
detection, and the crosses show SDSS ugriz data of the bright counterpart. The stars, squares, and triangles are the 90, 150, and 180 
fim. ISOPHOT data, respectively - filled symbols mark detections, whereas unfilled symbols show an upper limit. The best fit GRASIL 
SED is overplotted as the red solid line, and the type with the redshift if indicated. 



turbed morphology and NGP03 is the clearest case of an 
interacting system. Assuming typical K-band mass-to-li^ht 
ratios of star- forming spiral galaxies (|Bell De JondboOll : 
iGil de Paz et al] I2OOOI ) the estimated stellar masses range 
from 4 to 15 X 10^° i-^- close to a mature galaxy 



(jCole et al.ll200lh . The exception in this class is NGP24, a 
very nearby, 34 Mpc, star- forming dwarf galaxy with an ab- 
solute brightness of a mere Mr = —19.6. 

Should this set of galaxies, ignoring NGP24 from now 
on, be characterized as normal quiescent spirals or star- 
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Figure 3. Optical R -\- I combined images of the faint sample of the ISOPHOT NGP source fields, where multiple and/or faint and 
distant counterparts are identified. The large circle of one arcmin radius shows the location of the ISO detected source. If one of the 
counterpart candidates is a disturbed or interacting galaxy at z < 0.3, it is shown in the inset, and its x^y position in the main image is 
given. 



bursts? The SFRs derived by both GRASIL and ERR03 
fits are in the range of 7 — 26 Mq yr~^, i.e. higher 
than ordinary local spirals by factors of few, but still lower 
than typical strong nuclear starbursts. The ERR03 mod- 
els, Table [2l show them to be cirrus dominated, but with 
significant contributions from warmer dust in star-bursts. 
Determination of the specific star- format ion rate (SSFR) 
also place our galaxies somewhere in between quiescent 
galaxies and starbursts: SSFRs are \og{SFR/Msteiiar) ^ 
—9.8, si milar to local Ha selected star-forming g alaxies 
(see e.g. iGil de Paz et al.1 l2000l : ISaiina et al] l2006h . The 
ISOPHOT sample SSFR values are clearly higher than 
those of local galaxies in general, even the 'blue' ones: 
most of the stellar mass in the local Universe peaks at 



\og{SFR/Msteiiar) —11.6 wMlc the blue cloud is a t 
around \og{SFR/Msteiiar) ^ -10.2 fe.g. lKauffm"aiinll2004l ). 
Note some ambiguity in the term "starburst" : some studies 
find "starbursts" while others find galaxies with "slightly 
elevated star- formation" compared to normal local spirals, 
while actually meaning similar kind of galaxies. 

The peaks of the rest-frame SEDs are in the range 
100 to 150 [im. According to the ERR03 modelling the 
galaxies have cold dust temperatures in the range 15 
to 20 K - these values refer to the dominant dust 
population, the fits are not fitted with a single dust 
species and temperature. Dust masses are quite uni- 
formly - 10^ Mq. IRAS log(/^(100iJ,m)//^(60M,m)) and 
Spitzer log(/iy(160^Lm)//iy(70^Lm)) colours of the unambigu- 
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Figure 3 - continued 



ous sources, as calculated from the best-fit models are in 
the range -0.3 to -0.5, and -0.3 to -0.7, respectively. These 
values are typical to the infrared properties of the more qui- 
escent FIR-selected galaxies, rather than active ones, in the 
nearby U niverse (e.g. Dale et al. 2001; Sanders et al. 2003|; 
bale et al.ll2005l) . Hence, the bright ISOPHOT galaxies in 
this field appear to be gas rich galaxies with clearly elevated 
star- formation at about the LIRG-class limit, but without 
evidence in our data for strong (nuclear) star-bursts. From 
optical data alone, e.g. with their quite uniform g — r^ 0.75, 
our sample is indistinguishable from other local relatively 
high mass "blue-cloud" galaxies. As they are close to the 
red- sequence, it could be speculated that they are moving 
towards it, the effects of which can be seen in their IR and 
morphological properties. 

The bright FIR selected sources in this survey and oth- 
ers tend to be (nearby) cold galaxies, whereas the local 



ULIRGs typically have warmer SEDs. However, at higher 
redshifts where LIRGs and ULIRGs start to dominate the 
energy budget of the Universe, Spitzer FIR selected galax- 
ies in fact app ear to have colder SEDs feheng et al.l l2QQ7l : 
ISvmeonidis et al. 2009, e.g.), more similar to the majority 
of local FIR-bright galaxies rather than local ULIRGs. This 
makes detailed studies of the physical conditions in local FIR 
selected samples of IR-bright non-ULIRGs very motivating. 



4.2 The nature of confused counterparts 

Among the rest of our FIR targets, the faint sample, there 
are several cases where there either are optional counter- 
parts, or strong reasons to suspect a blend of two or more 
bright galaxies. The (likely) properties of this set can be 
summarized as follows: 

A blend of two or more spiral galaxies at redshifts 
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Figure 4. SEDs of galaxies in ISOPHOT fields containing more than one good candidate for the FIR counterpart. For clarity, we overplot 
only Sc, M82 and Arp 220 SEDs normalized to the 150 fim fiux at fixed redshifts. See text for discussion on individual cases. 
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Table 3. The 'confused sample' of the ISOPHOT counterparts. We classify the sample according to the likeliest scenario for confusion: 
type-1 - sum of several bright galaxies; type-2 - sum of bright and faint galaxies; type-3 - sum of several faint galaxies; type-4 - one of 
several possible counterparts. 
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z < 0.3 is responsible for the ISOPHOT detection in 9 source is possibly blended with a lower-z counterpart, up to 
cases (NPG02, NGP05, NPG07, NGPll, NGP12, NGP14, a maximum of 36 per cent of our NGP targets include an 
NGP18, NGP23, NGP25). The photometric redshifts and IR-galaxy at ^ > 0.4. 
GRASIL SED types of the individual counterparts are es- 
sentially identical to the unambiguous cases in Sections 13.21 
and 14. 1[ though the optical colours show a slightly larger 
range than the unambiguous bright sample. Since we are 
not able to distribute the FIR flux to blended counterparts 
with spectroscopic redshifts, we will not attempt to model 
their physical characteristics in detail. However, we did run 
the ERR03 models on several cases where there e.g. were two 
spiral galaxies with reasonably secure photometric redshifts: 
the cirrus and starburst luminosities and dust characteristics 
and star- format ion rates again come out to be very similar 
to the bright unambiguous set of FIR sources. They are of- 
ten at slightly higher redshifts, but due to dividing the FIR 
fluxes to two or more sources, the Ljr and SFR values end 
up very similar. Thus, in total 14/22 or 64 per cent of the 
FIR sources in the ISOPHOT EBL project NGP fields are 
low-redshift z < 0.3 moderately star-forming galaxies. The 
discussion presented in Section [4.11 above applies directly to 
all these sources, the majority detected in our survey. 

Four other blended FIR sources (NGP06, NGP15, 
NGP16, NPG19) also include a contribution, possibly a ma- 
jor contribution, from a bright spiral or starburst. In these 
fields, however, we have also identified fainter red galax- 
ies which might well be higher redshift (U)LIRGs also con- 
tributing the FIR flux. 

4.3 Higher redshift counterparts 

There are four fields (NGP04, NGP09, NGPIO, NPG21), 
where the only likely counterpart appears to be a higher 
redshift, z > 0.4, (U)LIRG, or a sum of them. Higher red- 
shift counterparts thus account for a minimum of 18 per 
cent of the FIR sources. Including ones where a higher-z 



4.4 Comparing to previous surveys and models 

Overall we found 23 per cent of the NGP fields to have 
a uniquely identified target, while another 41 per cent are 
identified as blends of bright galaxies. In the various FIR- 
BACK and Lockman H ole surveys, typically ^ 50 per cent 



of sources are identified f'Pennef eld et al.l 20051: iTavlor et ahl 
I2OO5I : lOvabu et all 1200 5: Rodig hiero et al.1 l2005l ). and the 
rest remain unidentified. The fractions are similar to ours 
even though the methods of identification have been quite 
different: most of the above works considered only a bright 
sub-set, typically half, of their total source catalogs, and 
they all employed mid-IR and radio data in the identifica- 
tion which we have not done. 

If the general fractions of identified and unidentified 
sources are similar to other works, what about the charac- 
teristics of the counterparts? We identified at least 64 per 
cent of all targets as fairly normal cool or cold IR lumi- 
nous spiral galaxies at z < 0.3, regardless of confusion is- 
sues. If the cases of likely bright and faint galaxy blends 
are included, the fraction rises to 82 per cent. Similarly, the 
fraction of contributions from higher redshift (U)LIRG type 
galaxies is in between 18 and 36 per cent. It is important 
to note that these cases are not different in their FIR prop- 
erties from the rest of the NGP ISOPHOT sources. Again, 
the results are very consistent with other surveys. In general, 
previous ISOPHOT FIR follow-up surveys have found a bi- 
modal distribution of a large, approximately three-quarters 
majority of low-redshift quiescent galaxies an d a minority 



highe r redshift ULIRGs at z = 0.4 — 0.9 (e.g. 'Sa iina et al 



2003: Patris et al.ll2003l: iTav 



120051 : iDennefeld et al.|[2QQ5l : 



or et al.ll2005: Rodighiero et al 



Saiina et al.ll2006l ). In the stud- 
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ies probing somew hat fainter flux levels of 100 mJy (e.g. 
lOvabu et al]l2005l ) than ours and FIRBACK's -150 mJy, 
the fraction of sources at redshift higher than 0.3 starts 
to ris e slightly. T he d ust te mper atures we find ar e simi- 
lar tolSaiina et all (e.g. bood ) and lPatris et al.1 (|2003l ) while 
iTavlor et al.l (|2005l ) derive slightly higher temperatures in 
the range 20-40 K. Early results from the FIR surveys with 
AKARI also suggest large contributions from local z < 0.1 
galaxies (Matsuura et al. 2007. : .Malek et al . 2009). 

Spitzer MIPS results also result in a broadly equiva- 
lent picture: the deeper MIPS sources at 70 |J,m (down to 4 
mJy) appear to be starb ursts at z=0.1-1.2 with a mean and 
median at z — 0.5 (e.g. ISvmeonidis et al.ll2008h . while the 
160 |J,m sources (to 100 mJy) are mostly "starburst galax- 
ies" ai z ^ 0.2 though generally not of the active, warm 
SEP type of M82 an d ARP 220, but with cooler SEDs (e.g. 
iFraver et aL 2006b; Wen et al.|[2007h . 

We do identify a difference to previous follow-up sur- 
veys, however. Four out of five of our bright unambiguous 
counterparts, as well as more than half of the identifiable 
blended counterpart options are clearly disturbed or inter- 
acting galaxies. In contrast, Sajina et al. (e.g. 2006) found 
that their counterparts to ELAIS Nl survey do no t show 
evidence of interactions, and also ,iOvabu et all (|2005h cite a 
small minority as interacting sources in the Lockman Hole 
survey. We are not sure what the reason for this difference 
is, unless it is just that our optical data is somewhat deeper 
than theirs making morphological classification easier. 

The majority of the galaxy counterparts in "deep" FIR 
surveys discussed in this paper are in fact very similar to 
the large sample of much brighter galaxies at > 2Jy de- 
tected as part of the ISOPHO T 170 ^tm Serendipity Survey 
(jStickel Klaas k Lemk3l2007f ). 

4.5 Models 

Previous FIR surveys have been described t o a r easonable 
accuracy by models of e.g . Lagache et al.l (|2004l l and re- 
cently by 'Rowan- Ro binso] 



20091 ). The models have been 



specifically constructed to fit the source counts at vari- 
ous wavelengths, and also the level and spectral shape of 
CIRB. They are broadly consistent with the bimodal dis- 
tribution of FIR sources described above. For example the 
160 ^im Spitzer counts are dominated by "cirrus" galaxies 
brighter than 80 mJy, and by M82 -type starbursts fainter 
than this limit according to the Rowan- RobinsonI (|2009l l 
models, and at 150 mJy, more or less at our survey limit, 
quiescent galaxies should outnumber starbursts by a factor 
of three, consistently with the nature of ISO sources found 
in the NGP fields. However, a recent FIR/sub-mm follow- 
up of BLAST at the longer wavelengths of 250 and 500 ^im 
(jPve et al.|[2009l ) find that these same bimodal models are 
a poor fit to the th eir data, as are also early AKARI counts 
of iMatsuura et al.l (|2007|), highlighting the fact that the ex- 
act nature of galaxies emitting their peak at 100 to 200 ^tm, 
where the CIRB also peaks, remains of great importance. 



4.6 Confusion at FIR wavelenghts 

In our NGP- field follow-up most of the FIR sources, 17/22 or 
approximately 70 per cent, cannot be identified unambigu- 
ously with a single optical counterpart, though there clearly 



are different classes of ambiguity. Of these 17 ISOPHOT de- 
tections 9 are cases where the most likely counterpart of the 
FIR detection is a sum of two or more fairly bright {K < 16 
mag) and nearby {z < 0.3) star-forming galaxies, and in an- 
other four cases one such bright nearby spiral appears to be 
a component in the confused FIR source. Further four fields 
do not have any nearby galaxies in them, leaving the expla- 
nation for the FIR flux to either be a single higher redshift 
ULIRG type source, or a sum of several of these - there is 
no way to differentiate between these cases with the present 
data. 

Thus, based on the present data, about 60 per cent of 
whole sample are definitely blended objects, with another 20 

S ier ce nt possibly so. Similarly to our survev. lPennefeld et al.l 
2005|l find 50 per cent of their identified sources multiple, 
and 'Saiin a et al.] (|2006l l note that up to 50 per cent of the 
FIRBACK sources are not ideally matched with a single 
counterpart. 



4-6.1 Effect on source counts 

How does the blending and confusion affect interpretations 
of FIR populations? In addition to making identifications 
and SED fits of individual galaxies difficult and sometimes 
suspect, it also affects the FIR source counts themselves. 
As an example. Fig. [5] shows the raw 180 iiim source counts 
from our ISOPHOT EBL project (jjuvela. Mattila &: Lenik3 
'2000^) as well a s raw FIRBACK/ELAIS counts from 
Dole e t al.l (j200l[ ). The solid symbols are original points, 
whereas the open symbols show the effect of blending. In 
the left panel we plot the NGP counts, with the following 
modifications: green open squares are a case where the 9 ob- 
jects definitely consisting of bright blends are divided into 
two discrete equal- flux sources; blue diamonds show a case 
where each detection in the faint sample is divided into two 
discrete sources of equal flux; and magenta triangles depict 
a case where each detection in the faint sample is divided 
into two sources where one has 2/3 and the other 1 /3 of the 
flux. 

To see the effect with somewhat larger number statis- 
tics, the second panel shows the counts from the whole 
ISOPHOT EBL survey. In this case we do not have the 
information on possible blending of each individual source, 
so we make a statistical estimate. The blue open diamonds, 
equivalently to the left panel, show the average of 10000 re- 
alizations of the original counts, where randomly selected 60 
per cent of the FIR sources are replaced with two sources of 
half of the original flux each. The error bars reflect the stan- 
dard deviation of the simulated counts. The triangles show 
the case where the random 60 per cent are divided into a 2/3 
and a 1/3 flux source, as above. We remind the reader that 
there is no difference in the flux distribution of the blended 
and unique sources, justifying the selection of split sources 
from the whole population. Since FIRBACK/ELAIS surveys 
saw very similar fractions of non-unique counterparts (e.g. 
iDennefeld et al. I I2OO5I : ISaiina et al. II2OO6II, we ran the same 
test on the FIRBACK 175 |xm catalog, and the result is 
shown in the right panel of Fig. O where the symbols and 
colours are equivalent to the first two panels. In addition, 
we ove r plot t he raw MIPS 160 ^im counts of iFraver et al I 
(|2006bl . I2OO9I ) with green circles. Only those points in the 
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FIRBACK and MIPS counts which have completeness of 85 
per cent and above are plotted. 

Though the faintest plotted raw count bins are affected 
by incompleteness, it still is clear from Fig.[5]that the reshuf- 
fling of flux typically takes some 30 to 60 per cent of the 
abundant 100 to 400 mJy sources and shifts them into flux 
bins in between 50 and 200 mJy. The resulting differential 
counts move down broadly by a factor of ^ 1.5, and steepen 
the faint end in the presented examples where 60 per cent 
of sources are affected . We ran the statistical blending cor- 
rection with a variety of affected percentage of sources, and 
a range in the fractions and number of underlying sources 
corresponding to a FIR detection: e.g. when 30 per cent of 
sources are affected in a way that 1/5 of a source flux is given 
to a new faint source, the necessary downward correction of 
counts is still approximately a factor of 1.2. Where splitting 
is more even the correction is larger, and rises to factors of 
few if three or more underlying sources are concerned. Only 
if the affected source population is below 20 per cent, are 
the required corrections negligible compared to the errors of 
the observed counts. 

Interestingly, the MIPS counts fall closer to the 
blending-corrected ISOPHOT counts and are not compat- 
ible with the uncorrected ones. The explanation is in the 
confusion levels. There are two kinds of "confusion": one 
due to fluctuations of all the unresolved sources below the 
detection Umit and the other due to the effective beam size 
not resolving properly sources too close to each other. One 
can talk about confusion "noise" in the first case, and e.g. 
"blending" in the second case, and of the respective confu- 
sion limits as a "photometric confusion limit" and a "source 
density criterion" often given in terms of beams per source 
in the latter case. Both effects are non-trivial as they obvi- 
ously depend on the shape of the source count slope both 
abo ve and below th e confusion limits and other factor s (se e 
e.g. [Vaisanen et al.ll200ll : IPole et al.ll2003l : iTakeuchi Ishiil 
l200r and references therein). For the cases here, the MIPS 
160 |J,m beam of ~ 40" is smaller than the respective C200 
ISOPHOT beam of ^ 90^^ size, and assuming a Euclidean 
slope this means that the confusion "noise" is expected to 
be 3 times higher for ISOPHOT, and eve n more for s teeper 
cou nts. Indeed, f or a s et of source counts iDole et al.1 (|2003l ) 
and ljeong et al.l (|2006l l estimate the confusion "noise" limits 
to be in the region of 50 and 160 mJy for MIPS 160 ^tm and 
ISOPHOT 180 ^Lm, respectively. These values correspond 
roughly to ~ 40 or more beams per source, which is a limi t 
one should use with steep counts (|Vaisanen et al.1 l200lh . 
Hence, the ISOPHOT counts we are considering clearly 
reach right to the lev el of exp ected confusion, whereas only 
the two deepest Fraver et al J ([2009) points at < 100 mJy in 
the Fig. [5] approach the MIPS confusion limit. 

An added complication in confusion estimates is clus- 
tering of sources which increases the confusion limits with 
an increasing severity as the source counts get steeper - 
ITakeuchi Ishiil i 2004l l show how the effect can be factors 
of a few or more for super-Euclidean slopes. An attempt to 
estimate the effect quantitatively is beyond the scope of this 
paper, but as we have clearly seen, we often detect close 
pairs or groups of likely counterparts within the FIR beam 
size, so the confusion limits cited above are likely lower lim- 
its. It thus should come as no surprise that the ISOPHOT 
FIR galaxy counts, both ours and others in the literature. 



have been significantly affected by confusion, as probably 
are the very deepest MIPS FIR counts. One should also 
be aware that any upward bumps in source count slopes 
when approaching the c onfusion reginie, are likely due to 
confusion noise effects (i Vaisanen et al.ll200lL e.g.) and not 
real sources, especially in the case of super-Euclidean counts 
such as these. 

In summary, the effects of confusion and blending on 
source counts, their likely too high normalization and differ- 
ences in slopes, must be properly taken into account when 
modelling the count slopes with galaxy populations. Lumi- 
nosity functions based on FIR counts will be affected as well. 



5 SUMMARY 

We have presented follow-up observations in the NGP fields 
of the ISOPHOT EBL Project. This FIR survey comple- 
mented the other major ISO FIR projects of ELAIS, FIR- 
BACK and Lockman Hole, by detecting sources in three 
bands, and also b y being able to determine the absolute 
level of the CIRB (jjuvela et al .|[2009). 

The NGP fields consist of 1.64 sq.deg area and 25 unique 
FIR sources in the 90 |J.m band and 22 in the 150 and 180 
\im. We imaged these fields with the NOT in the optical 
and NIR, and employed SED fitting techniques as well as 
morphological analysis to determine the counterparts of the 
FIR sources in the l' radius ISOPHOT error circles. 

Only five sources were securely identified with a single 
bright nearby galaxy. One of these is a local dwarf galaxy, 
but all the others are IR-luminous Ljr ^ 10^^ , fairly massive 
0.3-l.Om^ star- forming spirals at redshifts of z < 0.2. Their 
star- formation rates range from 7 to 26 Mq yr~^, dust tem- 
peratures are 20 K and dust masses Md ^ 10^ ^0- 
Such cold and fairly local galaxies have been found before 
in ISO and Spitzer FIR surveys selected in the 100-200 
[im range. We note in addition that in our data most of 
these appear to have disturbed morphologies showing signs 
of present or past interactions . 

Nine more FIR targets were securely identified with 
multiple galaxies. It turns out, however, that these cases 
have essentially identical physical characteristics with the 
first single-galaxy counterpart group, with the slight dif- 
ference that their redshift range extends to z ^ 0.3. Thus 
2/3 of our FIR galaxies are definitely relatively normal IR- 
luminous cold star- forming galaxies at z < 0.3. 

Half of the remaining cases, nearly 20 per cent of the 
total, also have a bright star- forming galaxy in the field, but 
to be explained the FIR flux needs additional contribution 
from an optically fainter galaxy. Finally, the second half of 
the remaining cases have only fainter galaxies in the area. In 
all these cases we were able to find two or three good candi- 
dates for LIRGs or ULIRGS in the range z = 0.4-0.8, which 
could explain the FIR flux. Without further spectroscopy 
the exact identification remains ambiguous however. 

Given the large number of blends, we tested the of- 
ten ignored aspect of how the confused sources affect the 
FIR source counts themselves. By using both our results 
and those from previous surveys, we showed that significant 
reshuffling of FIR fluxes is very likely to happen: a large 
fraction of apparent FIR sources at the 200-500 mJy range 
are in fact bound to be double or triple sources each with 
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Figure 5. Uncorrected differential source counts in the 180 \xm band, normalized with a Euclidean slope, are shown with solid symbols 
for the NGP fields at left, all ISOPHOT EBL Project areas in the middle, and ELAIS Nl and N2 fields on the right. Open symbols 
show effects of blended sources on the source counts; see text for details of confusion and blending effects. A model pre diction from 
iLaga che et al] (|2004l ) is shown for reference as the red curve. The right panel also shows the Spitzer/MIPS 160 ia.m counts of lFraver et al.l 
(|2006bl . I2009t) . which reach a factor of 2-3 fainter than the ISOPHOT counts displayed, as green circles. 



fluxes in the range 50-300 mJy. This has the effect of both 
decreasing the normalization of the FIR counts by a factor 
of ^ 1.5 — 2 and also steepening the count slope somewhat, 
especially at the faint end. This has to be taken into account 
when modelling FIR source counts. 

Finally, many of the NGP ISOPHOT FIR sources ap- 
pear to be part of significant galaxy concentrations, pairs, 
groups or perhaps even clusters. What this tells about the 
FIR sources and the reasons why they are IR-luminous, will 
be investigated in more detail in another paper. 
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APPENDIX A: DISCUSSION OF INDIVIDUAL 
ISOPHOT FIELDS 

Al Bright unambiguous galaxies 

NGPOl The bright (i?=16.59, i^=13.77) counterpart is 
identified as 2MASSXJ13412738+4042166 with a redshift of 
z = 0.088. The optical data show a disk dominated galaxy 
with wide arms, typical to late spirals, bar structure iden- 
tified with GALFIT and a somewhat disturbed morphol- 
ogy at the southern edge, typical to e.g. merger remnants. 
HYPERZ fits perfectly a fully evolved (nominally 15 Gyr 
old) GRASIL Sc template, the ERR03 models fit a cirrus- 
dominated SED, while the SDSS spectrum shows strong Ha 
emission. The FIR points alone - detections at 150 and 170 
/im, but only an upper limit at 90 /im - suggest a strong cool 
cirrus component. We also note that there is another bright 
star- forming galaxy at exactly the same redshift just outside 
the FIR detection circle, 2.5 arcmin (or 250 kpc projected 
distance) away from the adopted counterpart. 

NGP03 The bright (i^=17.12, K=14.20) counterpart 
is identified as 2MASSXJ13422216+4022017 at z = 0.131. 
It has two or three fainter galaxies clearly in close inter- 
action. GALFIT modelling shows a large bulge-dominated 
disk of 6.9 kpc scale length with wide spiral arms and the 
subtraction of bulge and disk components reveals a ring- 
Uke structure extending towards two of the companions, and 
possibly bar-like structure at the nucleus. Again, the SDSS 
spectrum shows strong star formation. The FIR SED is now 
warmer (detected at 90 /xm but not at 180 /im) and our best 
fitting SED template is, indeed, that of a prototype inter- 
acting starburst galaxy NGC 6090, further corroborating its 
starburst nature. Note that this template has more power 
in the cooler FIR range than an M82 SED would have. The 
ERR03 model is similar to NGPOl. 

NGP20 The bright (i?=17.75, i^^=14.80) galaxy in the 
detection circle has no SDSS spectra, but it does have a 
strong HYPERZ fit (x^ 0.2 at z = 0.137) with a fully 
evolved Sc template, consistent with the ERR03 fit as well. It 
has an undisturbed disk and a small but bright bulge in the 
optical. There is another bright galaxy of early type just out- 
side the field 1.1 arcmin away (2MASX J13494050+3907555) 
which has a known redshift of z = 0.143 making the photo- 
metric redshift of the first galaxy, the adopted counterpart, 
very plausible. 

NGP22 The only bright galaxy (i?=17.70, 
i^=14.52) in the vicinity is the edge-on disk galaxy 
SDSSJ135054.71+385847.2 at a redshift of z = 0.086. The 
disk is somewhat warped, though there are no obvious 
major partners nearby. The SDSS spectrum shows Ha, 
though not much else, and the overall SED fits reasonably 
the evolved Sc template at that redshift, though the opti- 
cal/NIR SED would be better fit with an M82 template. 
This galaxy is 2-3 times less massive than the previous 
four objects, and the SFR predicted by both GRASIL and 
ERR03 is also smaller, but still a little higher than ordinary 
quiescent galaxies. ERR03 again fits a cirrus dominated 
galaxy. 

NGP24 The counterpart galaxy of this FIR source is 
the most nearby galaxy in our sample, UGC 08793 at a 
redshift oi z — 0.0081 and it extends over 1.5^ on the sky. Its 
IR luminosity is more than 2 orders of magnitude less than 
the other bright sample galaxies above. The SDSS spectrum 
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of the center shows a typical late type spiral spectrum, NED 
classifies it as Sd, and our overall SED shape is well fit with 
a template of local Scd dwarf galaxy NGC 6946. This is the 
only case where ERR03 suggests a pure cirrus spectrum. 
It has numerous bright HII regions, several small satellite 
galaxies, or giant HII regions, and the disk appears slightly 
warped. 

A2 Confused fields 

NGP02: There is a significant concentration of K ^ 15—17 
mag galaxies within and just outside the ISOPHOT area - 
there are four galaxies with K < 16 mag and a further four 
with K < 17 mag. None have spectroscopy available and 
none of the galaxies within the detection circle produce satis- 
factory SED fits when ISOPHOT data is included. However, 
the four brightest galaxies all have optical/NIR photomet- 
ric redshifts of z ~ 0.28, along with further three fainter 
ones. They include a tidally disrupted interacting pair, two 
Sb-type SEDs, as well as two bright early type galaxies. A 
combination of two or three of the disk galaxies at this red- 
shift is well able to produce the required FIR flux. An exact 
identification of NGP02 is thus ambiguous, though a mixed 
origin is very likely. 

NGP04: There are no bright {K < 16, R < 20) galaxies 
in this field. The FIR detection is 130^^ away from an X-ray 
detected galaxy cluster RX J1342.8+4028 at z = 0.699. The 
HYPERZ fits to SDSS and RIK data find a large concentra- 
tion of z ~ 0.6 galaxies in the ISOPHOT area, which most 
likely belong to the cluster. Specifically, there are three red 
galaxies within the detection circle which fit quite well an 
Arp 220 SED at those redshifts - and the fits are even better 
if the FIR is a sum of two of them. Two of these sources, 
including the brightest object in the field, are clearly dis- 
turbed objects. A sum of faint galaxies is thus the likeliest 
counterpart for the ISOPHOT detection. 

NGP05: The brightest galaxy is a morphologically dis- 
turbed barred late type spiral; there are no spectra avail- 
able. The galaxy fits reasonably an evolved Sc template at 
z ~ 0.25, though it would need ^ 2 times more FIR flux 
to fit well the FIR points. Other galaxies may thus con- 
tribute to the FIR flux; and indeed there are 2 other galaxies 
(K < 16.6 mag) with photometric redshifts at z ~ 0.2, best- 
fit SEDs (when fit is done excluding the ISOPHOT points) 
of early and late type spirals and disk morphologies. A con- 
fused counterpart is thus likely. 

NGP06: The brightest galaxy is 2MASSX 
J13430669+4014314 with an SDSS redshift of z = 0.163. 
Both the optical spectrum and optical/NIR SED indicate 
an Sb-type, but the SED including all three FIR points 
does not fit well any of our templates at that redshift - the 
NGP 6090 template gives a reasonable fit at the 2a error 
levels of the FIR fluxes, but the optical spectrum is not 
that of a starburst. It is possible that the FIR detection 
come from some combination of a warm (contributing to 
the 90 /im band) and a cold object (the 180 /xm flux). There 
are also two fainter interacting pairs in the detection circle 
which both would reasonably fit an Arp 220 template at 
z ^ 0.5 — 0.6, contributing to the 90 and 150 /xm fluxes in 
that case. Th e identification thus remains amb iguous. 

NGP07: lJuvela. Mattila Lemk^ (l2000l ) list this ob- 
ject with two positions in the 90 /im catalogue (NGP 07 



and 08, separated by 20 arcsec) and a single detection in 
150 and 180 /xm, suggesting a complex origin for the FIR 
flux. There are 4 fairly bright galaxies in the area, though 
only one K < 16 mag. By far the brightest appears to be an 
early-type spiral or lenticular while the three others are also 
all early type disks. None of these alone fit the SED tem- 
plates over all the wavelength range. A confused FIR source, 
in between z ^ 0.25 — 0.40 based on photometric redshifts, 
is highly probable. 

NGP09: The two brightest galaxies in the region, only 
one is i^^ < 16 mag, are disk galaxies, but we cannot fit 
well either of them and the FIR flux is too strong to be 
associated even with their sum. There are, however, several 
very red galaxies giving acceptable fits to ULIRG templates 
in a wide redshift range of z ^ 0.2 — 0.8. We also note that 
this is one of three sources in the sample of 22 NGP FIR 
targets which do not have quality flag values of ^ ^ 3 in any 
FIR band. 

NGPIO: There are no galaxies brighter than K — 16 
mag in the field and no object provides an acceptable HY- 
PERZ fit. There are, however, half a dozen galaxies redder 
than R — K > 4 in the area including two EROs, which 
could contribute in combinations at z ~ 0.5 — 1.0 where 
their optical/NIR photometric redshifts place them. 

NGPll: No single galaxy fits well the template SEDs, 
though there are two K ^ 15.5 disk galaxies in the field, 
the sum of which can explain the FIR points if both are Sc- 
types at their photometric redshift ol z ^ 0.3. They belong 
to an apparent concentration of galaxies partially inside the 
ISOPHOT area. We find, however, another fainter red dis- 
turbed galaxy which could contribute as a higher redshift 
ULIRG. 

NGP12: Four galaxies, all K ^ 16 mag or brighter, 
produce excellent exponential profiles with GALFIT and 
good optical/NIR SED fits to spiral or starburst templates 
in the range z = 0.10 — 0.19. One, or even two, of them ap- 
pear to be interacting with another galaxy. However, none 
of them fit our full SED templates alone, though a sum of 
them could explain the required FIR flux. To complicate 
matters, there are four fainter red R — K ^ 4 galaxies in 
the field, any of which fit reasonably well a full ULIRG 
template in the range z = 0.4 — 0.55. The very ambigu- 
ous aspect of this source is hi ghlighted by the fact that 
Juvela, Mattil a Lemke (2000) list this target as two ob- 
jects (NGP 12 and 13) separated by 23^^ in the 90 /xm cata- 
logue, while there is only one detection in the 150 /xm list. 

NGP 14 There are two bright nearby spirals in the 
area. The brighter one (i?=17.32, K=14.19) is 2MASSX 
J13473443+3931515, with an SDSS redshift of z = 0.16. 
The optical spectrum looks to be of a typical star-forming 
galaxy, and the GRASIL interacting starburst NGC 6090 
template gives a perfect match to the overall SED. The op- 
tical image reveals a disk galaxy with a large asymmetric 
arm. This case would fall into the bright unambiguous class 
were it not for the second spiral (K— 15.22) which also fits 
very well our templates, an evolved Sc-galaxy at a photo- 
metric z — 0.08. Both galaxies must contribute to the FIR 
flux. 

NGP15: This is one of the most ambiguous cases. 
There is a tight group of 5-7 galaxies (all fainter than 
K ^ 15.5 mag), including two interacting pairs, 25^^ from 
the FIR location. Two or three of these galaxies together 
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would fit reasonably the FIR points with starburst templates 
at z ^ 0.3. Additionally, however, there are three very red 
sources elsewhere inside the FIR area. If their elliptical-like 
optical/NIR colours are combined with the FIR flux, they 
fit well ULIRG templates at z ^ 0.6 — 0.9. And, in fact, these 
sources belong to a distribution of the largest concentration 
of EROs found in this follow-up survey, mostly lying just 
outside the ISO detection circle. NGP15 is one of the three 
sources with lower quality flag values of ^ < 3. 

NGP16: There are two separate 90 /xm detections just 
22" from each other (NGP16 and NGP17), that are con- 
nected with a single detection at longer wavelengths. A sig- 
nificant concentration of galaxies is seen in the area. There 
are two disk galaxies in the center of the area, with early type 
spirals as best-fit SEDs at photometric redshifts of z ^ 0.2, 
but even a sum of these would not yet explain the FIR fluxes. 
Intriguingly, there are also three individual fainter and red 
(R — K ^ 4) galaxies in the detection area which each give 
a good fit to starburst or Arp 220 SEDs at redshifts ranging 
z — 0.1 — 0.5. One of these is a clear interacting galaxy, and 
another a FIRST radio source whose radio flux is consistent 
with the fitted Arp 220 SED. Any one of these could be the 
true counterpart, unless it is not the confused sum of the 
brighter galaxies, or a combination of all. This is a highly 
ambiguous case. 

NGP18: This field is within a galaxy cluster ZwCl 
1346.9+3931. Three bright disturbed spirals are found in 
the ISOPHOT area, and several more just outside, all 
with spectroscopic SDSS redshifts of z = 0.13 — 0.14. 
The brightest (i?=17.02, K=14.19) of the spirals, 2MASX 
J13490845+3917219, is a star forming galaxy with Seyfert- 
like emission lines, and the overall SED actually fits well 
an evolved Sc template. On the other hand, since the other 
spirals in the area must also contribute to the FIR flux (opti- 
cal/NIR SEDs are fit with normal spiral and even starburst 
SEDs), we chose to include NGP18 in the confused sam- 
ple, rather than the bright unambiguous sample. It is clear 
though that the FIR counterpart (s) lies at the redshift of 
z^O.13. 

NGP19: The ISOPHOT circle does not have very 
bright galaxies, though there are many at both sides of 
K ^ 16 mag. The 90 /xm flux is the strongest of the whole 
sample, making this our only source with a clearly declin- 
ing SED longwards of 100 /xm. Morphologically, according 
to GALFIT, two of the brighter ones are face-on disks with 
stellar-PSF and/or strong bulge components. Both fit well 
an Arp 220 template at z ^ 0.35, though without spectra 
it is impossible to speculate more about starbursts or ob- 
scured AGN, or to decide between them. Contributions from 
other galaxies cannot be ruled out either: the second bright- 
est galaxy in the field is a warped interacting disk galaxy, 
it does not fit any template by itself, and there are three 
fainter very red galaxies in the area giving reasonable fits to 
ULIRG or ERO templates in redshift ranges of z ~ 0.3 — 0.6. 

NGP21: There are no galaxies brighter than K — 16 
(the bright object seen in Fig. [3] is a star), and none of the 
sources give good SED fits. The brightest source is a disk 
galaxy, but no template fits the FIR points. There are also 
three very faint EROs in the field which could conceivably 
contribute if they have ULIRG type SEDs, though none of 
them alone. This FIR source is the third ISO target with a 
lower quality flag value of ^ < 3. 



NGP23: This field again holds a significant concentra- 
tion of bright galaxies. The brightest (i?=16.73, ^^=13.60) 
galaxy has the GALFIT profile of an elliptical or strongly 
bulge dominated disk galaxy and it fits reasonably an overall 
Sb-type SED. The two next brightest galaxies have SDSS 
spectra {z ~ 0.20) showing star formation, and they are 
disks, as are two other bright galaxies, all with optical/NIR 
spectral shapes consistent with spirals. Most of these galax- 
ies must contribute toward the total FIR flux, and the iden- 
tification thus remains ambiguous. 

NGP25: This case is somewhat ambiguous though 
there is a single bright galaxy (i?= 16.65, 13.61) in the 
ISOPHOT area: it appears to be an early type galaxy (Sa 
template fits the optical/NIR SED), with some evidence of 
structure beneath a GALFIT de Vaucoulers profile. It does 
not have a spectrum, but there is a galaxy of the same size 
and brightness 47^^away to the NE just outside the detection 
circle. Its SDSS redshift is z = 0.118 and it has a clear disk- 
like GALFIT profile, and the best fit SED is an Sb-type. 
Contributions from (at least) both of these galaxies would 
be needed, however, to explain the FIR flux. An Sc-type 
SED could fit the overall SED alone but the morphology of 
neither galaxy favours this option. 



Table Al. Photometry of the unambiguous bright counterparts from our /, K-band NOT imaging, our ISOPHOT observations (Juvela 2000), while the ugriz data are from the 
SDSS archive. Spectroscopic SDSS redshifts are given with 3 decimal places, whereas our own Hyperz-photometric redshifts with two decimal places. Note that the RA and DEC refer 
to the optical/NIR position of the FIR counterpart, not the FIR detection. 
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i 
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90|am 


150 |nm 


180 |nm 


redshift 


NGPOl 


13h41m27.4s 


40°42'16'' 


18.80±0.07 


17.46ib0.02 


16.80ib0.03 


16.59zb0.01 


16.41zb0.02 


16.03=b0.01 


16.24zb0.10 


13.74zb0.01 


< 170 


490 zb 60 


500 zb 60 


0.088 


NGP03 


13h42m22.2s 


40°22'02'' 


19.44dz0.09 


18.11ib0.01 


17.34dz0.01 


17.10zb0.01 


16.93±0.01 


16.47zb0.01 


16.69zb0.03 


14.17±0.01 


170dz50 


200 di 70 


< 300 


0.131 


NGP20 


13h49m34.9s 


39° or 30'' 


20.09ib0.14 


18.70ib0.02 


17.93dz0.01 


17.75±0.01 


17.53±0.02 


17.13dz0.01 


17.28di0.02 


14.77±0.01 


< 150 


150 ± 60 


200 zb 80 


0.14 


NGP22 


13h50m54.7s 


38°58'4r' 


19.84ib0.18 


18.60ib0.04 


17.84dz0.02 


17.70±0.01 


17.45±0.02 


17.06±0.01 


17.25±0.08 


14.49ib0.01 


< 120 


320 ± 110 


420 zb 110 


0.086 


NGP24 


13h52m35.0s 


38°42'18'' 








18.11di0.01 




18.04dz0.02 




15.32dz0.05 


< 130 


260 zb 60 


280 zb 80 


0.008 
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Ob 



O 
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Table A2. A sample of photometry from all the NOT data in the NGP fields as described in Sections 2.2 and 3.1 of the Paper is shown here. The full catalog of 15714 sources is 
available electronically from CDS {http://cdsarc.u-strasbg.fr/) - contact the first author for access before the paper appears in print. Note that the catalog contains sources also from 
outside the ISOPHOT error circles. The columns list the ALFOSC field id, a running object number N within that field, RA and DEC, then the SExtractor CLASS (CL) galaxy/star 
classification parameter from ALFOSC data, I and K SExtractor total magnitudes (e.g. Rt), 2.1" aperture magnitudes (e.g. R2.7), and their errors (e.g. eRt), the SDSS DR7 ugriz 
magnitudes and their errors, and finally the IP indicates whether the object is within 60^^ of a given numbered ISOPHOT FIR NGP source. 



R2.7 



K2.7 



6^2.7 



347 
348 
349 
350 
351 
352 
353 
354 
355 
356 
357 
358 
359 
360 
361 
362 
363 
364 
365 
366 
367 
368 
369 
370 
371 
372 
373 
374 
375 
376 
377 
378 
379 
380 
381 
382 
383 
384 
385 
386 
387 
388 
389 
390 



205.340103 
205.340302 
205.340500 
205.340698 
205.340805 
205.340805 
205.341095 
205.341797 
205.341995 
205.341995 
205.342407 
205.342499 
205.342896 
205.342896 
205.343002 
205.343307 
205.343399 
205.343399 
205.343597 
205.343800 
205.343903 
205.343994 
205.343994 
205.344498 
205.344604 
205.344894 
205.345093 
205.345200 
205.345200 
205.345500 
205.345596 
205.345596 
205.346405 
205.346893 
205.346893 
205.347107 
205.347198 
205.347600 
205.347900 
205.348100 
205.348206 
205.348297 
205.348404 
205.348404 



40.733978 
40.722130 
40.694592 
40.700760 
40.696388 
40.739101 
40.697529 
40.660759 
40.667881 
40.683739 
40.657299 
40.680012 
40.730228 
40.728432 
40.744381 
40.661411 
40.666248 
40.688320 
40.724831 
40.712950 
40.697868 
40.654282 
40.729069 
40.696041 
40.659061 
40.692242 
40.683819 
40.673759 
40.721111 
40.702600 
40.648869 
40.651150 
40.705742 
40.705860 
40.708519 
40.726608 
40.717892 
40.716520 
40.708080 
40.721030 
40.703239 
40.675339 
40.679119 
40.710121 



0.77 
0.37 
0.41 
0.75 
0.03 
0.01 
0.16 
0.65 
0.03 
0.01 
0.44 
0.51 
0.04 
0.38 
0.55 
0.03 
0.36 
0.42 
0.04 
0.00 
0.03 
0.38 
0.44 
0.94 
0.27 
0.44 
0.64 
0.47 
0.16 
0.00 
0.31 
0.37 
0.66 
0.62 
0.75 
0.03 
0.98 
0.00 
0.39 
0.00 
0.47 
0.46 
0.03 
0.16 



22.84 
24.95 
23.33 
23.50 
19.47 
20.31 
22.15 
22.48 
19.70 
22.11 
23.54 
24.56 
19.24 
24.33 
21.68 
18.21 
23.06 
24.91 
22.02 

0.00 
19.24 
24.00 
23.64 
22.70 
22.01 
24.71 
23.72 
23.43 
21.87 

0.00 
22.76 
24.71 
23.41 
23.55 
22.72 
19.41 
17.64 

0.00 
24.24 

0.00 
22.35 
24.26 
17.57 
22.30 



22.85 
24.87 
23.35 
23.50 
19.96 
21.51 
22.76 
23.20 
20.12 
22.43 
23.68 
24.47 
19.66 
24.19 
23.20 
18.77 
23.12 
24.69 
22.20 

0.00 
19.67 
24.07 
23.91 
22.65 
22.21 
24.65 
23.78 
23.67 
22.12 

0.00 
23.12 
24.63 
23.47 
23.66 
22.75 
20.39 
17.75 

0.00 
24.25 

0.00 
22.84 
24.30 
18.73 
22.49 



23.09 
25.69 
22.26 
22.15 
19.09 
19.87 
21.47 
22.56 
18.96 
21.25 
23.40 
23.58 
18.55 
23.70 
21.30 
17.55 
22.42 
24.61 
21.06 

0.00 
18.61 
23.81 
23.41 
22.16 
21.59 
23.91 
22.23 
23.13 
21.54 

0.00 
22.11 
24.64 
23.29 
22.56 
22.35 
18.82 
16.05 

0.00 
23.57 

0.00 
21.55 
23.55 
17.12 
21.39 



22.93 
24.67 
22.30 
22.11 
19.55 
21.08 
21.92 
23.18 
19.32 
21.51 
23.71 
23.67 
18.91 
23.53 
22.50 



21.23 

0.00 
18.98 
23.79 
23.62 
22.11 
21.67 
24.06 
22.22 
23.19 
21.70 

0.00 
22.31 
24.66 
23.09 
22.59 
22.19 
19.65 
16.13 

0.00 
23.70 

0.00 
21.99 
23.50 
18.18 
21.61 



0.06 
0.22 
0.10 
0.10 
0.01 
0.03 
0.07 
0.07 
0.01 
0.05 
0.13 
0.21 
0.01 
0.23 
0.06 
0.00 
0.08 
0.28 
0.04 
0.00 
0.01 
0.22 
0.15 
0.05 
0.05 
0.23 
0.13 
0.14 
0.04 
0.00 
0.11 
0.23 
0.12 
0.16 
0.07 
0.01 
0.00 
0.00 
0.21 
0.00 
0.08 
0.23 
0.00 
0.07 



0.06 
0.36 
0.09 
0.10 
0.01 
0.02 
0.05 
0.07 
0.01 
0.04 
0.13 
0.25 
0.00 
0.20 
0.09 
0.00 
0.07 
0.30 
0.03 
0.00 
0.00 
0.19 
0.15 
0.05 
0.03 
0.29 
0.13 
0.12 
0.03 
0.00 
0.08 
0.31 
0.10 
0.12 
0.05 
0.01 
0.00 
0.00 
0.21 
0.00 
0.06 
0.21 
0.00 
0.04 



0.13 
0.74 
0.07 
0.05 
0.01 
0.03 
0.06 
0.13 
0.01 
0.04 
0.20 
0.15 
0.00 
0.22 
0.08 
0.00 
0.08 
0.37 
0.03 
0.00 
0.00 
0.31 
0.21 
0.05 
0.06 
0.19 
0.06 
0.19 
0.05 
0.00 
0.10 
0.38 
0.18 
0.11 
0.09 
0.01 
0.00 
0.00 
0.19 
0.00 
0.06 
0.21 
0.00 
0.05 



0.11 
0.53 
0.06 
0.05 
0.01 
0.02 
0.04 
0.11 
0.00 
0.03 
0.23 
0.21 
0.00 
0.19 
0.08 
0.00 
0.07 
0.40 
0.02 
0.00 
0.00 
0.25 
0.20 
0.05 
0.03 
0.30 
0.05 
0.13 
0.04 
0.00 
0.06 
0.56 
0.12 
0.08 
0.05 
0.01 
0.00 
0.00 
0.22 
0.00 
0.05 
0.18 
0.00 
0.03 



0.00 
0.00 
0.00 
18.80 
17.32 
0.00 
19.18 
0.00 
16.01 
18.64 
0.00 
0.00 
15.69 
0.00 
0.00 
0.00 
0.00 
0.00 
17.90 
19.48 
16.10 
0.00 
0.00 
0.00 
0.00 
0.00 
18.70 
0.00 
0.00 
19.47 
0.00 
0.00 
0.00 
0.00 
19.09 
16.34 
13.60 
18.90 
0.00 
19.67 
0.00 
0.00 
15.15 
19.31 



0.00 
0.00 
0.00 
18.91 
17.64 
0.00 
19.33 
0.00 
16.25 
18.98 
0.00 
0.00 
15.99 
0.00 
0.00 
0.00 
0.00 
0.00 
17.86 
19.68 
16.28 
0.00 
0.00 
0.00 
0.00 
0.00 
18.83 
0.00 
0.00 
19.69 
0.00 
0.00 
0.00 
0.00 
19.06 
16.80 
13.63 
18.96 
0.00 
19.47 
0.00 
0.00 
16.04 
19.20 



0.00 
0.00 
0.00 
0.13 
0.09 
0.00 
0.26 
0.00 
0.03 
0.22 
0.00 
0.00 
0.02 
0.00 
0.00 
0.00 
0.00 
0.00 
0.07 
0.34 
0.02 
0.00 
0.00 
0.00 
0.00 
0.00 
0.20 
0.00 
0.00 
0.20 
0.00 
0.00 
0.00 
0.00 
0.21 
0.04 
0.00 
0.20 
0.00 
0.19 
0.00 
0.00 
0.02 
0.14 



0.00 
0.00 
0.00 
0.19 
0.06 
0.00 
0.28 
0.00 
0.02 
0.20 
0.00 
0.00 
0.01 
0.00 
0.00 
0.00 
0.00 
0.00 
0.07 
0.37 
0.02 
0.00 
0.00 
0.00 
0.00 
0.00 
0.17 
0.00 
0.00 
0.39 
0.00 
0.00 
0.00 
0.00 
0.21 
0.03 
0.00 
0.19 
0.00 
0.30 
0.00 
0.00 
0.01 
0.24 



23.98 
0.00 
0.00 
0.00 
20.49 
24.04 
24.13 
0.00 
21.70 
24.14 
0.00 
0.00 
26.06 
0.00 
22.35 
20.78 
0.00 
0.00 
22.69 
0.00 
21.16 
0.00 
0.00 
24.27 
25.38 
0.00 
0.00 
0.00 
24.17 
0.00 
0.00 
0.00 
0.00 
0.00 
24.02 
21.44 
21.55 
0.00 
0.00 
0.00 
0.00 
0.00 
19.43 
24.28 



24.20 
0.00 
0.00 
0.00 
20.01 
21.36 
23.66 
0.00 
21.18 
22.91 
0.00 
0.00 
20.85 
0.00 
20.94 
19.37 
0.00 
0.00 
22.98 
0.00 
20.28 
0.00 
0.00 
23.33 
23.29 
0.00 
0.00 
0.00 
22.64 
0.00 
0.00 
0.00 
0.00 
0.00 
23.88 
20.57 
19.40 
0.00 
0.00 
0.00 
0.00 
0.00 
18.25 
23.87 



22.62 
0.00 
0.00 
0.00 
19.55 
20.78 
23.37 
0.00 
19.97 
21.87 
0.00 
0.00 
19.41 
0.00 
20.60 
18.37 
0.00 
0.00 
22.34 
0.00 
19.38 
0.00 
0.00 
22.29 
22.23 
0.00 
0.00 
0.00 
22.09 
0.00 
0.00 
0.00 
0.00 
0.00 
22.61 
19.55 
17.99 
0.00 
0.00 
0.00 
0.00 
0.00 
17.73 
22.75 



22.26 
0.00 
0.00 
0.00 
19.66 
20.61 
22.14 
0.00 
19.40 
21.69 
0.00 
0.00 
18.93 
0.00 
20.87 
17.90 
0.00 
0.00 
21.42 
0.00 
19.01 
0.00 
0.00 
22.05 
21.95 
0.00 
0.00 
0.00 
22.11 
0.00 
0.00 
0.00 
0.00 
0.00 
21.97 
19.22 
16.59 
0.00 
0.00 
0.00 
0.00 
0.00 
17.42 
22.13 



21.39 
0.00 
0.00 
0.00 
19.66 
20.87 
22.53 
0.00 
19.23 
20.83 
0.00 
0.00 
18.43 
0.00 
21.06 
17.71 
0.00 
0.00 
21.30 
0.00 
18.71 
0.00 
0.00 
22.25 
21.39 
0.00 
0.00 
0.00 
22.05 
0.00 
0.00 
0.00 
0.00 
0.00 
21.97 
19.26 
15.81 
0.00 
0.00 
0.00 
0.00 
0.00 
17.26 
21.17 



0.93 
0.00 
0.00 
0.00 
0.15 
1.67 
0.98 
0.00 
0.44 
1.72 
0.00 
0.00 
1.61 
0.00 
1.58 
0.22 
0.00 
0.00 
0.40 
0.00 
0.27 
0.00 
0.00 
0.87 
0.99 
0.00 
0.00 
0.00 
0.93 
0.00 
0.00 
0.00 
0.00 
0.00 
0.92 
0.61 
0.20 
0.00 
0.00 
0.00 
0.00 
0.00 
0.10 
0.92 



0.40 
0.00 
0.00 
0.00 
0.03 
0.07 
0.28 
0.00 
0.11 
0.26 
0.00 
0.00 
0.07 
0.00 
0.15 
0.02 
0.00 
0.00 
0.17 
0.00 
0.04 
0.00 
0.00 
0.19 
0.28 
0.00 
0.00 
0.00 
0.11 
0.00 
0.00 
0.00 
0.00 
0.00 
0.32 
0.09 
0.01 
0.00 
0.00 
0.00 
0.00 
0.00 
0.01 
0.29 



0.17 
0.00 
0.00 
0.00 
0.03 
0.06 
0.31 
0.00 
0.06 
0.16 
0.00 
0.00 
0.03 
0.00 
0.17 
0.02 
0.00 
0.00 
0.15 
0.00 
0.03 
0.00 
0.00 
0.11 
0.15 
0.00 
0.00 
0.00 
0.10 
0.00 
0.00 
0.00 
0.00 
0.00 
0.16 
0.05 
0.01 
0.00 
0.00 
0.00 
0.00 
0.00 
0.01 
0.16 



0.20 
0.00 
0.00 
0.00 
0.07 
0.09 
0.19 
0.00 
0.05 
0.23 
0.00 
0.00 
0.03 
0.00 
0.38 
0.02 
0.00 
0.00 
0.11 
0.00 
0.03 
0.00 
0.00 
0.14 
0.18 
0.00 
0.00 
0.00 
0.17 
0.00 
0.00 
0.00 
0.00 
0.00 
0.15 
0.08 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.01 
0.16 



0.48 
0.00 
0.00 
0.00 
0.33 
0.59 
0.94 
0.00 
0.17 
0.56 
0.00 
0.00 
0.11 
0.00 
2.29 
0.05 
0.00 
0.00 
0.50 
0.00 
0.13 
0.00 
0.00 
0.72 
0.40 
0.00 
0.00 
0.00 
0.72 
0.00 
0.00 
0.00 
0.00 
0.00 
0.71 
0.40 
0.01 
0.00 
0.00 
0.00 
0.00 
0.00 
0.06 
0.36 



